A simple strong solution to non-linear stochastic HJB PDEs: an 
application to the portfolio model 



Abstract. We devise a simple and general method for solving non-linear 
stochastic Hamilton- Jacobi-Bellman partial differential equations. We apply 
our method to the portfolio model. 
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1. Introduction 

This paper overcomes major obstacles in the area of stochastics and 
mathematical finance. We provide simple strong solutions to the (non- 
linear stochastic) Hamilton- Jacob i-Bellman partial differential equation HJB 
PDE (and other stochastic PDEs). In so doing, we transform the stochastic 
HJB PDE into an ordinary differential equation ODE. It is well-known that 
stochastic HJBs are generally non-smooth and that the existing methods rely 
on the methods of viscosity and minimax weak solutions (see, for example, 
Crandall and Lions (1983)). 

We apply our method to the portfolio model (the stochastic factor model). 
It is well-known that a criticism of the portfolio model is that the optimal 
portfolio does not depend on wealth for an exponential utility function (the 
most common utility in finance). Consequently, in this paper, we resolve 
this issue. That is, we show that the optimal portfolio depends on wealth, 
regardless of the functional form of the utility. 

2. The method 

We consider this function V (x, y) ; it can be expressed as V (j3x, y) , where 
(5 is a shift parameter with an initial value equal to one (see Alghalith (2008)). 
Define g = (3x; differentiating V (g,y) with respect to j3 and x, respectively, 
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yields 



V p = V g x, 
V x = V g (3. 



Thus 



£ = Uv. = ™. (1) 

Vg x x 
The second derivatives of V (g,y) with respect to B and x, respectively, are 



VfiP = VggX 2 , 



Therefore 



Vxx X 2 Tr X 2 Vrr , \ 

— = -9 => V xx = — (2) 

Vhr B 2 B 2 V ' 



Similarly, we can rewrite V (x, y) as V (x, j3y) . Define / = By; differentiating 
V (x, f) with respect to B and y, respectively, yields 
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Vp = V f y, 



Thus 



(3) 

The second derivatives of V (x, f) with respect to (3 and y, respectively, are 



Vpp = Vffy 2 , 

Vyy = V ff P 2 . 



Therefore 



Vyy _ v P%1 u . 



Differentiating flTJ with respect to y yields 



my 



xy 



V, 



xV 3 



X 



Py 



xy 



(5) 



Differentiating ([3]) with respect to y yields 



Substituting (jSJ) into (jSJ) , we obtain 



3. Practical example: the portfolio model 

We provide a brief description of the portfolio model (see, for example, 
Detemple (2012), Alghalith (2009) and Castaneda-Leyva and Hernandez- 
Hernandez (2006), among others). Thus, we have a two-dimensional Brow- 
nian motion I^W^Wj^J , J- S j on the probability space (fi, J 7 ,, P) , 
where {J r s } t<s<T is the augmentation of filtration. The risk- free asset price 

T 

/ r(Y s )ds 

process is 5*0 = 6' , where r (Y s ) e Cj> (R) is the rate of return and Y s 
is the economic factor. 

The risky asset price process is given by 



dS s = S t {fi(Y s )dt + a (Y s ) dW l s } , (8) 
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where \i (Y t ) and a (Y t ) are the rate of return and the volatility, respectively. 
The economic factor process is given by 



dY s = b(Y s )dt + pdW] + v 7 ! - P 2 dWP ,Y t = y, (9) 

where \p\ < 1 is the correlation factor between the two Brownian motions 
and b (Y s ) e C 1 (R) with bounded derivative. 
The wealth process is given by 



X* = x + J {r (Y s ) X* + (/i (Y s ) - r (Y s ) vr s )} ds + J n s a (Y s ) dW l s , (10) 



where x is the initial wealth, {-n t , J~ s } t< <T is the portfolio process is the 

T 

consumption process, with E f n^ds < oo . The trading strategy n s G A (x, y) 

t 

is admissible. 



The investor's objective is to maximize the expected utility of terminal 
wealth 



V(t,x,y) = SupE[u(X^\J t }, (11) 
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where V (.) is the value function, u (.), bounded and strictly concave utility 
function. 



The value function satisfies the Hamiltonian-Jacobi-Bellman PDE 



V t + r (y) xV x + b (y) V y + -V m + 



Sup^nfa 2 (y) V xx + [n t (// (y) - r (y))) V x + pa (y) n t V xy j = 0, 

V(T,x,y) = u(x). (12) 

Thus we obtain the following well-known HJB PDE (suppressing the nota- 
tions) 



V t + rxV x + (/i-r) ix* t V x + l -o 2 vfV xx + bV y + hf m + pair* t V xy = 0, (13) 



where the asterisk denotes the optimal value and we suppressed the notation. 
Substituting (}!]) , , , l|g)) and (|Tj) into (jlSj) yields 



K+^Va+lM-HTTj +-tr vr t — — + +- — ^+pcx7r t p - 

x 2 /r y 2 y z xy 
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Using the above procedure we can easily show that V t = (3Vp/t and thus the 
above equation becomes 

f ^ + r M - r) + W + = o. 

t x 2 p V 2 y z xy 

The solution to the above equation is equivalent to solving an ODE. Since 
V (.) is differentiable with respect to the shift parameter j3, the solution is 
classical (strong), even if V (.) is not smooth in x, y or t. 
Setting j3 at its initial value, we obtain this boundary condition 

v -i + rv e+{ , - r)*t% + iov^ + + i% + ^ei+M = o. 

t x 2 y 2 y z xy 

Thus the optimal portfolio is given by 

n * = (£ - r) Vp _ pa(V^ + Vfs) 
1 cr 2 x 3 V/3f} crx 3 yVf3p 

We note that, in contrast to the previous literature, that the optimal portfolio 
depends on wealth, regardless of the functional form of the utility, even for 
an exponential utility function. Clearly, this is a more realistic result. 
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